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Abstract 



A complete description of the linearized gravitational field on a flat background is given 
in terms of gauge-independent quasilocal quantities. This is an extension of the results from 
^ , . Asymptotic spherical quasilocal parameterization of the Weyl field and its relation with 

' Einstein equations is presented. The field equations are equivalent to the wave equation. A 

I generalization for Schwarzschild background is developed and the axial part of gravitational 

■ field is fully analyzed. In the case of axial degree of freedom for linearized gravitational field 

the corresponding generalization of the dAlembert operator is a Regge- Wheeler equation. 
Finally, the asymptotics at null infinity is investigated and strong peeling property for axial 
waves is proved. 

^ ■ 

CT'. 1 Introduction 

{3J}, We show that seemingly complicated linearized Einstein equations on a Schwarzschild background 

can be analyzed in terms of gauge invariants. The obtained invariants decouple, in a natural 
way, into axial and polar parts keeping symmetry with respect to the interchange of the null 
' coordinates u and v. The invariant y describing axial degrees of freedom, corresponding to 35'o, 

fulfills Regge- Wheeler equation because axial part of the corresponding component of the Weyl 
field is gauge-invariant^ On the other hand the polar part the Weyl field is not gauge-invariant. 
However, all components of the Weyl field may be "corrected" in such a way that we obtain 
invariants which substitute linearized Ncwman-Penrose scalars. 

In we have shown how the gauge-invariant quantities x, y describing unconstrained degrees 
of freedom of the gravitational field arise in a canonical formalism. Here, we concentrate on their 
relations with linearized Weyl tensor. In the case of a flat background we present, in Theorem 1, an 
explicit relation between linearized Weyl tensor and the invariants. We continue this analysis for 
the case of a Schwarzschild background and we show an analogous relation but only for the axial 
degree of freedom described by y. Finally, we examine this result in view of the so-called peeling 
property and we show that axial part of the linearized gravitational field obeys strong peeling. 

This paper is organized as follows: In the next Section, some preliminary notions and results 
for the flat background are introduced. Section 3 contains a generalization for the Schwarzschild 
background, in particular, gauge invariants and Einstein equations for the axial part of the gravi- 
tational field are presented. In Section 4 we discuss the relation between invariants and linearized 

*Partially supported by a grant KBN Nr 2 P03A 047 15 and CNRS Orleans. E-mail; Jacek.Jezierski@fuw.edu.pl 
^More precisely, some components in axial part of the Weyl field are gauge invariant but not all of them and 
they have to be "corrected" by metric terms (see Section W) to become gauge independent. 
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Riemann tensor. Section 5 is devoted to the investigation of the asymptoties at nuU infinity for 
the solutions of Regge- Wheeler eqiiation and its connection with peeling property. To clarify the 
exposition some of the technical results and proofs have been shifted to the appendix. 

2 Description in null coordinates for flat background 

We present in this section some standard results about linearized gravitational field with nontrivial 
extensions not only in a notation but also in the framework. 

2.1 Minkowski metric in null coordinates 

Let us consider the flat Minkowski metric of the following form in spherical coordinates 

rj^^dy^dy" = -di^ + dr^ + r'^{de^ + sin^ ed(f>'^) . (2.1) 

The Minkowski space M has a natural structure of a spherical foliation around null infinity, more 
precisely, the neighbourhood of looks like S'^ x M2. We shall use several coordinates on M2: 
t, r, p, V, u. They are defined as follows 

p := u := t — r v := t + r . 

Let us fix the null coordinates (u, v) together with the index a corresponding to them. 

The coordinates on a sphere we denote {x^), {A = 1,2), {x^ = 0,x^ = (p) and the round metric 

on a unit sphere by 7ab ( ^ Asdx'^dx^ = d^^ + sin^ Odcp'^). Let us also denote by A the laplacian 

o 

corresponding to the metric 7 ab- Moreover, we use the symbol "||" for the covariant derivative 

on S'^ with respect to the induced metric tjab- 

For convenience we need also some more denotations: p = = ^71:;^, p,a = P^Sa where 
£u := h, Eti := r]°'^SaSb = 1- We define s" := ij^^Sb and one can check that = 1, = — 1, 

The explicit formulae for the components of Minkowski metric can be denoted as follows 

-2 ° 1 I I 

VAB = P IaB, 'nab = -^\Eab\, r?aA = 

where E^u = = E^y and = 1 = —Eyu and 

ry^^daj^d.x'" = r]abdx"-dx'' + r]ABdx'^dx^ = -dudv + p-^{d0^ + sin^ Odcj)^) . 
Similarily, the inverse metric has the following components 

^ p2 °ab ^ ^ab ^ _2|£;«i'| ^ r?"^ = 
where = = E"" and E""" = 1 = -E^". We shall also need the derivatives 

V^^,a = 2/0£a77^^ , TjAB ,a = -"^P^aVAB 

and finally the nonvanishing Christoffel symbols are the following 

T"aB = pe°-riAB , r^aB = -pSaS^B , T^BC 

where T'^bc are Christoffel symbols for the spherical covariant derivative "||" on S"^. 
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2.2 Riemann tensor in null coordinates 

The linearized Riemann tensor R^^xs defined in an obvious way in terms of the second derivatives 
of the linearized metric hf^u by the formula 

has the following components in null coordinates: 

'^Rabcd — had,hc ^bd^ac ^bcad ^acbd 

2RabcD — haD,bc — hbDMC + hbc,aD ~ hac,bD + 

+P£b {haD.c + hcD,a — haco) — P^a {hbD.c + ^cD,fc " hhc^o) 

2RAbCd = ^dylliab + f^bC\\A,d ~ Kd\\AC ~ hAC,bd + 

+peb {hdA\\C - hdC\\A - hAC\d) + p£d {hbC\\A - hbA\\C - hAC,b) + 

+P'nAC£^ {hbdM — had,b — hab.d) - '^P^Sb^dhAC 

2RABCd — hdAWBC + hBC\\A,d — hBd\\AC ~ ^AC||B,d + 2p£d (/lsC||A ^ /iylC||s) + 
+PVBCS"' {haA,d - hdA,a + had,A + 2pedhaA) - PVAC^'" {haB,d - hdB,a + had,B + 2pedhaB) 
2RabCD = haD\\C,b — ^b_D||C,a + ^6C||D,a " ^aC||D,b + 
+ 2peb {haD\\C — haCWo) + 2pea {hbC\\D — ^fcD||c) 
2RaBCD = /iylZ5||BC + hscWAD — hBD\\AC — ^AC||BD + 

+P11AC£'^ {hBDM - haB\\D - haD\\B) + PVBDS'" {hAC.a - ^aC||A - ^aA||c) + 
— PVBCS"' {hAD,a — haA\\D ~ /iaD||yl) " PVADS"' {hBC,a — haB\\C — ^aC||s) + 

(hBDVAC + h-AcVBD - h-ADVBC " hBcVAo) + 2p'^ e'' hab ivAcTlBD - VbcVad) 
We show in the sequel how the above formulae can be generalized for the Schwarzschild background. 

2.3 Ricci tensor in null coordinates 

The linearized Ricci tensor i?^^ := rj^'^Rs^Xv takes the following form in our coordinates 

2Rab = h^b.ac + ha'^ ^cb — hab'^c ~ h'^c.ab + haAM^^^ + hbA.a^^ — hab^^^ A ~ H ^ab + 
+peaH^b + pSbH^a + 2pe'' {hab,c - hacM - hbc^a) 

2RaB = h'^B.ab — haB''^ c + ha'^ ,cB — h'^c,aB + /ia'^||BA " haB^^^A + XB'^\\A,a ~ -;^H\\B,a + 



-psa (2h''B.b - h\^B) - 2pe''hbBM - 2p^eae''k 



bB 

C 



2Rab = {h^AWB + h^BWA)^^^ - h^aWAB - XAS'^a " 2p£"xAS,a + XA \\CB + XB \\CA + 
1 



-XAB^^^C + VAB 



^{HW^c + H'\) + 2pe-{H,a - ha^HA) + p\2s-e'Kb - H) 



where H := rj^^hAB and xab '■= ^ab ~ ^VabH. Some components of the linearized Ricci are 
also derived for Schwarzschild background in Appendix 



2.4 Gauge in null coordinates 

The gauge transformation 

splits in the following way 

hab > hab + Ca,6 + ^6,0 

haA > haA + £,a,A + £,A,a + 2pea£,A (2.2) 

hAB > hAB + ^A\\B + S,B\\A - 2prjAB£''£,a ■ 

There are also some useful formulae 

XAB > XAB + £,A\\B + S,B\\A ^ VAB^"^ \\C 

which are straightforward consequences of the previous ones. 

2.5 Invariants and vacuum Einstein equations 

Let us introduce the following gauge invariant quantitie^ 

ya := {k+2)haA\\B£^'' -{p'\A''\\CB£''''),a (2.3) 

y := 2p-\h,B\\Ae^''),aE'^'' (2.4) 
X := AH + p-W,a-i/ + 2e"eX6-2p-i£»/i,^ll^ (2.5) 

A(A+2)/i„f,-(A+2)[(p-2/i/||^),, + (p-2/i/||^)^a] 

+ [p-'{p-\^''uB)A,b+[p-'^p-'^^''\\AB)A^a (2-6) 

where e^^ is the Levi-Civita skew-symmetric tensor on a sphere {u — const., v — const.} such 
that p-^sinfe^^ = 1. 

The axial invariants are not independent they are related as follows: 

2p-2y^^^^ab ^ ( A + 2)y 



which is a simple consequence of the definition (2.4) 



If we assume that vacuum Einstein equations R^^i, — are fulfilled, we obtain the following 
equations for our invariants: 
The axial part reads as 



(p-V'^),a - 2p-^ Ra^'iibdS^'' = (2.7) 



^We leave to the reader an exercise to check that those quantities are gauge invariant (using (2.2)), however, in 
the Appendix N we show this property in a more general case. 
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or takes another form in terms of the quantity y"^ 

[p-\yafi - Yb.a)]'" + A + 2)y, = -2p-\ A + 2)RaB\\De'''' 
The polar part takes the foUowing form 

x'^^a^, - p2 A ( A + 2)x = 4(p-4 ^^^n^^).", + A ( A + 2) {'q'^^'Rab - v^^'Rab) = 

77'^^x,b = 4p-4^^^IIB^ = (2.9) 
- 2(p-2x) + Tj,,{p-^xy', = . (2.10) 



The left-hand side of the last equation (2.10) depends on all seven polar components of the Ricci 

tensor Rab, R'^^\\b, rj^^RAB-, R^^\\ba- Assuming that all of them are vanishing one can show 
that (2.10) is true. From eq. (2.1C) and (2.8) we conclude that the invariants Xab and y^ depend 
locally on x, y. More precisely, 

Xab = 2(p"^x)^ab - VabiP^^^yc 
y- = -2i?'^^(p-2y),6p' 

and the primary data (x, y) fulfills usual wave equation. 

(p"V)'\ + pAy = (p-ix)^^ + pAx = 
We describe in the sequel how the full Riemann tensor can be reconstructed from the invariants x, y. 

2.6 Quasi-local relations between gauge invariants and linearized Rie- 
mann or Weyl tensor 

It is convenient to use skew-symmetric tensor e""^ instead of density E°'^. It can be defined as 
follows 



v/|detr,,b|£™ = 1 



^1 detijabl 



-1 e' 



ab 



2E' 



ab 



One can show that the linearized Riemann tensor has the following (2 + 2) "spherical decomposi- 
tion" . In terms of our invariants it decouples into axial part 
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-pes HabCD — y 

RabcDWEs'^'e'''' = p'{p-'y)^, 
e^''RAB''d\\c^e/pHp-'y),b 

4p"* R^bCdWAoe^" = (p"Vrf),6 + (P"'yfc),d 
and polar part 



P V V tiABCD ^2 ^ABCD — X 

— Zp e e Habcd — P rj V ^abcd — —p V V tCaBcD ~ X 



(2.11) 
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P e ^ABCd\\E£ 



CE 



4p~* R^b'^dWAc = 

P~^{k+2)Rabc^\\D£ 



Xbd 

^P {P^ac),b£ 



(2.12) 



2p"2( A+2)i?^f,Ad = XM + (p ^x,d)^b + (p ^x,6),d 

In above formulae, as in the whole paper, we use extensively some operators on a unit sphere 
which become isomorphisms when we assume that mono-dipole part of the field vanishes (see here 
Appendix ^ and also |^] or Appendix B in j|] ) . The above equations contain the full information 
on ten independent components of the Weyl tensor up to the mono-dipole part of the field^. 
Moreover, one can easily check the "peeling" property at starting from the invariants x, y 
as a primary data. More precisely, assuming the following expansion 

(2.13) 
(2.14) 



X = Xip + X2P^ + X3/9^ + X4p^ 



and the same form for y 

„2 



y = yip + y2p 

we have 



,1 • N 2 

X„ = -Xip+ (-Xl - X2)p ■ 



1 



x„ = -xip 



-X3/9 - X4/9 



x„„ = 2xip ^ + 2(X2 - Xi) + . . . 
Xut, = ^^3/0^ + 3x4p* + . . . 

We summarize below in the table the relation of our invariants with the Newman-Penrose 
scalars and the Christodoulou-Klainerman-Nicolo decomposition of the Weyl tensor: 



Price 


Weyl 


C-K-N 


polar 


axial 


N-P 


asymptotics 




o ^ 
W V vB 




2 

P ^vv 


p^{p^^yv),v 




P''(x3,y3) 


^'_2 


o ^ 
Wu uB 


_2 

P ^icKN 


2 


p^{p^^yu)^u 


*4 


p(xi,yi) 




P e vVabv 


P ^ PcKK 


p2(p-ix),„ 


pHp-'y).v 


*1 


p''(x2,y2) 




p-^e''''Wabu^ 


^ i— CKN 


p2(p-lx),„ 


pHp~'y),u 


*3 


P^(xi,yi) 






PCKN; CKN 


2 

p X 


p^y 


*2 


P^(xi,yi) 



where W^abB = W^abB — bW^ abD (cf- TS transformation in appendix) and f_,\s is the 
linearized Weyl tensor. 

One can easily check that among all ^''s only p^^'^Q^""" = = x + iy fuUfills usual 

wave equation. The mono-dipole-free parts of the invariants x, y correspond to the unconstrained 
degrees of freedom of the linearized gravitational field. Moreover, they play a natural role of the 
"positions" in the reduced initial data set on a Cauchy surface ([|ll[, [l^ ). 
Remark The Teukolsky equations |2^ for 5*0, ^1^4 on a Kerr background seem to be quite strange 
as primary equations because they are not deformations of the usual wave equation when we pass 
to the asymptotically flat region. We would like to stress that there exists a generalization for 
the notion of x and y on a Schwarzschild background and both invariants fulfill a deformed wave 
equation - Regge- Wheeler for y and Zerilli for x (see jl^). 

^ The mono-dipole part is disscussed in ]l2) and it corresponds to the charges related to the Poincare group. 
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Theorem 1 The linearized Riemann tensor for the vacuum Einstein equations depends quasilo- 
cally on the invariants (x, y) which contain the full information about the linearized gravitational 
field. Moreover, the invariants x and y fulfill usual wave equation. 

In other words any mono-dipole-free solution (x, y) of the wave equation gives a Weyl field: 
1 o 



Wabcd = 
WABcd 



-p XEab^cd 



1 



-P y^AB^cd 



Wa^cdWB = -7;£cd£\p^{P ^^),b 



WaBcd\\E£ 



BE 



:£cdP^{P V),a 



(2.15) 



o 1 4 1 4 

Wc d\\AB = -^P ^cd = -^P 



[P ^^),cd - -^Vcdip ^X)'\ 



j^r A ^BC 

W c BdllAC^ 



-p' [ec\p-'y)M + ed''ip-'y)M] 



and W p^xs fulfils field equations given by Bianchi identities. The formula e (|2.15| ) are also valid 
(outside origin) if we include mono-dipole part of the fields x and y (see |12[|). More precisely, 
a mono-dipole solution x — Amp + 12kp^, y = 12sp^ corresponds to the 10 Poincare charges (a 
monopole and three dipoles): m — mass, s — spin, k — center of mass and p — linear momentum 
which is related with center of mass by the relation p — (9„ + 9„)k. Moreover, the "charges" fulfill 

o o o o 

the following equations: dain = daS = daP = 0, Am = ( A + 2)p = ( A + 2)s = ( A + 2)k = 
which simply mean that m is a constant and p, s are constant dipoles. 

Remark One can check that the "exterior" bounds (10.19 -10.20) for the various null components 
of the Weyl field in ||l^ (p. R114-R115) are related to the most important radiative terms xi and 
yi which appear in a^^^^, P^^^^j Pckn and ctckn but the fall-off conditions on Uqkn and Pckn are 
slightly weaker. In particular it is not known how to improve the fall-off for Ockn in the case of 
curved space-time. This leads to the problem in the standard conformal approach where stronger 
fall-off on acKNn assumed. 



2.7 Hierarchy of asymptotic solution on scri for scalar wave equation 

Let us consider the wave equation in null coordinates (m, v) 

p-i(p-V)'"a+ A</5 = (2.16) 



and suppose we are looking for a solution of the wave equation (2.16) as a series (see 



if = ipip + ip2P^ + (PZP^ + 



(2.17) 



where each Lpn is a function on , dyipn = 0. 

If we put the series ( p. 17 ) into the wave equation (2.16), we obtain the following recursion 



9„V3„+i = - — [ A + (n - l)n\ipn ■ 
2n 



(2.18) 



The above formula is the same as equations 2, 3, 4 in [Q and equation 1.6 in ^ but written in 
a more elegant way. Let us finish this Section with the following Remark which is devoted to 

■^corresponding to VE'q'^ 
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so-called Newman- Penrose constants: 

o 

Remark . The kernel of the operator [ A + ^(^ + 1)] corresponds to the l-th spherical harmonics. 
The right-hand side of vanishes on the n — 1 spherical harmonics subspace. This means 

that the corresponding multipole in does not depend on u. In particular, for n = 3 we 

have quadrupole charge in the fourth order. The nonlinear counterpart of this object is called 
Newman-Penrose constant (see [ p^ , [Q, [T^). In particular, in our case, this constant is related to 
the quadrupole part of X4 and y4. 



3 Gravitational field on a Schwarzschild background 

In [mI it is shown how to generalize the gauge-invariant quantities x and y for the case of 
Schwarzschild background. We would like to investigate the following problems: 



1. What is the relation (similar to ( 2.15 )) between linearized Weyl tensor and the gauge invariants? 
and 

2. How this relation applied to the asymptotics for the solutions of deformed wave equation 
(Regge- Wheeler and Zerilli) clarifies a peeling property for the linearized Weyl tensor? 

In this section we summarize results of Adam Jankowski [To| who analyzed axial degree of 
freedom. We shall often use the same letters for the generalizations of the objects from flat to 
Schwarzschild background with the obvious identification when the mass m ~ Q. 



3.1 Schwarzschild metric in null coordinates 

Let us start with spherical coordinates r, t, 9, (p, hence the Schwarzschild metric has the following 
form: 

Tif^^dx^^dx" = - ^1 - — ^ dt^ + ^ \,„ dr^ + r"^ {dO^ + siii^ 9 dcj)^) (3.1) 

Spacetime M has a natural foliation with respect to spherical symmetry, more precisely, it splits 
into S'^ X M2, where 5^ is a two-dimensional sphere and M2 is described by coordinates r and t. 
This way spacetime (M^rj) with coordinates {u,v,9,(f)) splits into (^^,77152) and (-M2, 771^2)- we 
shall often denote by x"^ coordinates {9, (p) on a sphere, but on M2 we use null coordinates defined 
in terms of standard r and t as follows: 

u = t — r — 2mln(r — 2to) mG] — oo,+oo[ 

V = t + r + 2mln(r — 2m) ug] — oo,+oo[ (3-2) 

Radial curves ( 9 = const., <j) = const.) with fixed u = const, are null geodesies, similarly for 
V — const. 
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Our main interest will concentrate on the domain r > 2m marked by grey colour on the above 
Carter- Penrose diagram for Schwarzs child spacetime. 

Let us denote coordinates (m, v) by small Latin characters (a;°), spherical coordinates by capitol 

Latin characters {x^),{A = 1, 2), {x^ ~ 9,x^ = (p), and the unit round spherical metric by Iab- 

Iab dx'^dx^ = dO"^ + sin^ OdcfP . 

Let us also notice that M2 -L S*^, because rjaA — 0. 

Let us fix some more notation: symbol "||" denotes two-dimensional covariant derivative on S*^ 
compatible with induced metric rjAB , the equality 

VAB = Iab 

o o 

implies that the Christoffel symbols for rjAB and Iab are the same. Moreover, symbol A denotes 

Beltrami-Laplace operator for the unit metric Iab- We keep the same notation as in Section ^ 
with obvious generalization from flat background to Schwarzschild one. 

ri^i,dx^dx^ = riabdx°'dx'' + rjABdx^dx^ (3-3) 

rj^^dx'^dx" = " " v) '^"'^'^ '^^^'^^^ ^ '^"^^ ^^"^^^ 
From definition ( |3.2[ ) we have ^{dv — du) = k~^dr where by k we denote: 

fc:=l-^. 
r 

We shall analyze the region far away from the sources and assume that r > 2m which implies 
k > 0. Using ( |3.2D one can verify that 

r.a^-EaVk, (3.5) 
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where Sa is defined as follows: 

Ev = e„ = ^ki (3.6) 

£" = -k-i £" = fc-^ (3.7) 

and the indices are raised by the two-dimensional inverse metric rj"'^. This way we have introduced 
a unit-length radial vector 

Vab£ £ —V Sa£b — J- • 

One can also check that 

VAB,a = SaVAB (3.8) 

r 

and we have the following nonvanishing Christoffel symbols for the metric connection on M: 

r^Ba = -Si—ea, (3.9) 



\/k 

r%B=r?As— (3.10) 
r 

r^BC = ^ri'^'^{VDB,c + riDC,B - riCB,D) , (3-11) 

rri 1 

T\c = ^k^ {-S\sc-d\sb + s'^rjbc) ■ (3.12) 

Let us notice that T^bc are simultaneously Christoffel symbols for the induced two-dimensional 
metric r]\g2- More precisely, in usual angular coordinates {6, (j)) we have two nonvanishing compo- 
nents: = cot 9 and F^^^ = — sin 6*008 0. 

We can assign the symbols T^bc and T°-bc to covariant derivative on S"^ and M2 respectively, and 
denote (on the example of a covector) as follows 

Ub dh^a - ^fT^'ab (3.13) 

U\\B--=dBU-iFT^AB (3.14) 

Riemann tensor for the metric jj^uj^ equals to Weyl (Schwarzschild metric is a vacuum solution) and 
has the following nonvanishing components (up to the symmetries): 

C^BCD = ^{ScVBD - S^r]Bc) , (3.15) 

2?7? 

C\cd = -^{S^cVbd - S2vbc) , (3.16) 



m , 

Levi-Civita connection on S"^ has the following components of the two-dimensional Riemann tensor: 
- ' 1 

Similarly, the corresponding curvature for the metric r]ab has the form: 



C^aBb = -^Sir,ab. (3.17) 



R ^BCD = -^{S^CVBD - S^DVBC) ■ (3.18) 



2 2?T? 

R \cd = T{S\Vbd - S^Vbc) ■ (3.19) 
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Let us introduce Levi-Civita tensor for {S'^^7]\s2) and (Af2,^|Af2) respectively: 



tpAB rpab 



v/|det ricD\ ypet^feil 

where E"-^ are coordinates of tensor density ^ A ^ so the values are in set { — 1,0,1}. Similarly, 
E"^^ are coordinates for ^ A ^ , so we may write explicitly 

^^"^ = 4—n ' (3-20) 
sm d 

e"" = ^ . (3.21) 
k 

Let us notice that the metric "qab has a signature (+,—), hence its determinant is negative: 
det ?7q6 < . 

This is important when we raise indices in Eab'- 

Sab = VacVbds'"' = -\/\det ■qab\Eab = ^^kEab , (3.22) 
but on the corresponding sign is positive: 



SAB = VacVbds^" = ^/\AbUub\Eab ^ sinOEAB ■ (3.23) 

It would be useful to derive explicit formulae for first and second derivatives of e^'~^ and rj^'^ with 
respect to null coordinates on Al2 which are implied by (3.5): 

^AC^^^^^AC^^^ ^AC^^^^_^^AC^^ (3 24) 

e tafc = i-;j£aeb - —Vab)e'^^ , ?7 tab = i-^jSaSb - —Vab)v , (3.25) 



where obviously the objects e^'-^ and r/'^'-^ are scalars with respect to the covariant derivative (3.14). 



3.2 Gauge transformation for the linearized metric tensor h 

We shall analyze, from 2 + 2 decomposition point of view, the gauge transformation generated by 
infinitesimal diffeomorphism of AI for the linearized metric tensoi]^ 

h^.i, -> hf,^, + 2^(^.j,) . (3.26) 

Let us first split the covariant derivatives of the covector field — {£,a,£,A)- 

^a;b — db^a ~ ^S^^ ab — ^aXb 

Ca;A = (^A^a " ^F^^ aA = ia\\A + -— Ca^o 



iA;b = dbiA — ^fL bA — ^Atb H ^^^6 



r 



iA-B = dsiA - (.F^^ AB - £.f^^AB = £.A\\B —^f^'' 

^The standard linearization formulae have been shifted to the Appendix [a[ 
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and next apply them to gauge transformation of the tensor hf^j^: 

+ ^atb + ^bta , (3.27) 

2\/fc 

haA ^ haA + S,a\\A + iAta ^ —£aS,A-, (3.28) 

2^/k 

hAB hAB + iA\\B + ^B\\A —^eS^VAB ■ (3.29) 

Here the symbol o over two-dimensional tensor Iab denotes its traceless part: 

tAB-.^ tAB - ^VABV^^tcD ■ (3.30) 

In particular, we denote the traceless part of hj^^g by 
o 1 

XAB =hAB-= hAB - -jVABH , (3.31) 

where H :— rj^^hAB- 

From (3.29) we get the gauge transformation for xab'- 

XAB XAB + ^A\\B + ^B\\A " VAB^"^ \\c ■ 

Let us notice that gauge of xab depends only on S^a, hence it is not dependent on the part which 
is tangent to M2. 

The ten components of the tensor split naturally with respect to the 2 + 2-splitting of the 
spacetime M = x M2 into: 

• components hah in M2, 

• components Hab on S*^, 

• mixed components haA- 

However, for the description of the two degrees of freedom of the gravitational field one can divide 
ten components of the tensor h^^i, differently, into axial and polar part. They split as follows: 

• 7 polar components: hab, ha^\\A, X"^ bwacV^^ , H; 

• 3 axial components: haA\\B£^^ , X"^ b\\ac£^'^ ■ 

In this section we shall consider only axial part of the gravitational field. In particular, the gauge 
transformation for the axial components of the tensor h^^ reduces to: 

haB\\C£^^ — * haB\\C£^'~'' + (Cs| |C£"^'^)ta , (3.32) 

o 

B AC B AC , ^BC\ to oo\ 

XA ||Bc£ XA ||Bc£ + — :^ — (Csiice ) (3-33) 

(see also Appendix ^ . Let us notice that the gauge of all axial components depends only on 
so they are not dependent on the infinitesimal change of coordinates on M2. Moreover, we shall 
see in the sequel that using appropriate operators on 5^ one can produce from haB\\c£^'~'' and 
XA^ \\BC^^'~' a gauge invariant quantity. 
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3.3 Gauge invariants 

The following object 

Ya ■■= (A+2)(/i,^||c£^^)-(r\^^l|Bc£-^^)t,. (3.34) 

is invariant with respect to the gauge transformation described in the previous subsection. It is 
useful to introduce another object (cf. p4|): 

y := r^{hbA\\ce^'')tae'"'. (3.35) 



which is also gauge independent. Let us notice that definitions (3.34) and (3.35) are the same 



as formulae (2.3) and (2.4) respectively. This means that for axial and invariants there are no 
"background mass" corrections. This phenomena is no longer valid for polar invariant x (see p^). 
The explicit proof of gauge invariance property for and y is given in appendix The invariants 
Ya and y (introduced also in [ ^ ) are not independent but they fulfill the following identity 

rVatbe'" - (A +2)y (3.36) 

which is a straightforward consequence of the above definitions of the objects. The dipole part 
of y corresponds to the stationary solution of the field equations^ (sec [Q). We shall verify this 
property in the sequel analyzing asymptotics of the solutions at future null infinity . 

3.4 Field equations 

Let us consider linearized vacuum Einstein equations 

f^up • ^ p\U(7 ^" h u-.pa ^up ^ ^ ly'^i^P ^ ^ aC* uap) • (3.37) 

The axial part of linearized Einstein equations on a Schwarzschild background can be easily 
described in terms of invariants [ p^ (we present the details in appendix |^) 

£d"(rV)ta + rVd = -2rVsrf||z5e^^ = , (3.38) 

(rV)ta = 2r4?^^l|Bc£''^ = 0. (3.39) 

The above equations take the same form as in the case of the flat backgrounc|^. Moreover, they 
imply a second order hyperbolic equation for y: 

(r-2(r2y)*")ta + r-2(A +2)y-0. 
However, this equation is no longer a wave operator but Regge- Wheeler equation 

+ = (3.40) 

4 Gauge transformation of linearized Riemann tensor and 
axial invariants 

In one can find details of calculations which we would like to present in this section. Those 
calculations are related to the gauge transformations of r^^vpa '■— SR^^pa, 5Rp,upa and r^n := 5Rpu 
(see appendix]^. The linearized Ricci r^i^ is obviously gauge invariant but ri^i_yp„ and SR^^^pa- are 
in general gauge dependent. 

^We shall introduce them in the next subsection. 
^This phenomenon is valid only for axial part. 
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Usi ng th e formula ( 3.26 ) one can show that the gauge transformation of the hnearized Riemann 
tensor (A_/7) discussed in appendix ^ takes the foUowing form: 



'^^a-p^ vpa '^^a;p^ apu '^^a-.a^ pup ^^^-^-R Mcrp ' 

Let us denote the gauge term as 

5Rpupa 5R^i,p„ + gage((5i?pypcr) 
hence the corresponding components take the form: 

ga.ge{5RabED) = 0, 



(4.1) 



(4.2) 
(4.3) 



771 ^TTl 
ga,ge{dRABDd) = -^^miAVBD - ^IbVAd) + —{^BVDA - (.AVDB)td 

H 7 £d[^BVDA - KaVOB) , 



gage(^i?AbCd) = --^iVbd^AWc + Vdbic\\A) +VAcF{^)db, 



where 



db 



{-)tb + {-)td + ^{^Sb + ^bSd) - 2i^bd—e^£.a 
r r Lr r 



(4.4) 
(4.5) 



To construct axial part of the Riemann tensor SR^^pa we use some "spherical operators" and 
obtain the following gauge dependence: 



5RabeD\\A£^ ~* '''^ 5 RabeD\\A£^ + — {Vbe{S,D\\A£^ ) ,a ~ Tlae{S,D\\A£ 



.DA\ 



SRab d||ce SRab d\\c£ 



c ^AB 4,rn 



(4.6) 
(4.7) 



SRAbBd ^RAbBd —^{^A\\B + ^BWA — VAbS. \\D)Vbd ■ 



(4.8) 



SRa, 



bED 



5Ra 



bED ■ 



(4.9) 



Let us notice that the components of the linearized Riemann tensor (4.6-4.8) may be "corrected" 
in such a way that we obtain gauge-independent objects: 

T.2('A I OUR -DA^ah , 4m J, B ^AC\ _ 

r (A +'i)oHabcD\\A£ £ H e c\r xa \\bc^ )tb — 



Am 



r\A +2)r„fc,i5||^e^^£'^'' + ^eb^ir^XA^ \\BC^'''')tb + ^hoaUVb 



.da 



_ ^ AB ^TTl BCJ CJ AB 5?77< BC 

oRaB d\\c£ + -^"-dS||c£ —fAB d\\c£ + -^"ciB II C£ 
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° m o 2m 

oRAbBd H ^VbdXAB ^r^Bd H ^VbdXAB 

^RabED — 1'abED 

For completeness, we give here also the above above formulae in terms of the "true" linearized 
Riemann tensor r^^po- which is related to dR^^po- by formula (A.4). More precisely, the relation 
(A. 4) for various components of Riemann tensor gives the following: 

771 

TabED ~ SRabED , TAbED — SRAbED , fABEd — SRABEd ~ -^hAd^lBE 
TabeD — SRabeD ^hoaVbe , ^AbBd — SRAbBd ^hABVbd ■ 

One can easily verify that the attempt to exchange 6Rf_iuXp with r^^xp does not improve the 



invariance of ( |4.6| - 4.8). In particular, the "metric corrections" although in a different form will 
still occur. We prefer to use the tensor 5Rp,i,\p because it possesses all the symmetries of the usual 
full curvature tensor. 

The explicit formulae for the axial part of 5Rf^upa in terms of the metric h^v 

2SRabED — haDtb\\E + KEta\\D — haEtbWD ~ KotaWE + (4-10) 
H SbC^haDWE — 2/laBIID) ^ ^a{2hi,E\\D — 2/lfcD||e) 



-{2hBE\\A£d - 2hAE\\B£d + had\\A£°''rtBE - hadWB^'^VAE) 



2SRABEd — h-AdWBE — hsdWAE + ^BEWAtd — iT-AEWBtd - 

Vk, 
r 

+ — {hBdta£''r]AE - hAdtaS^riBE + hAatdllEBe" - hBa%driEA£°') + 

k 

+ 2--^e''edihAaVEB - hBaVEA) (4.11) 

'^SRabeD = haDtbe ~ hae^^^D ^ fT'bDtae + hbc4a\\D + (4-12) 
+ —^{haDte£b — hbOfe^a + hoeta^b — ^Dejb^a + hbe\\D£a — ^aeHD^b) 

2SRAbBd = hAdtb\\B — hbd\\AB — hAB^bd + hbB\\AXd + 

H — —{hAd\\B£b + hbB\\A£d — hbAWB^d — hBdWA^b) + 
^/k 

- — {hABtb£d + hABXdSb + hfdXbS^riAB + hbfXdS^VAB + 

+hbdXf£^lAB) + ^{hAB£b£d - hABVdb - 2hdbVAB) (4-13) 
enables one to express them in terms of the invariants ya and y: 

Ir's-'e^^SRabED = y (4.14) 
tHA +2)SRAB''diic£^'' + ^ir\A''\\Bce^%d - i£/( A + 2)(ry)t, - l^£/(rV)ta 

O 4:771 

= (A +2)yd + rVke'^iyaXd ~ ydXa) - —yd (4.15) 
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4r" 



AC 



(4.16) 



'fA ^DA^ab , bf2 B cAU\ 

[A +Z)dUabeD\\AS £ H [r XA \\BC'^ )tb 



AC\ 



[r(A+2)y]te + ^(r2y)te. (4.17) 



The above formulae describe the relation between axial part of the linearized Riemann tensor 
(corrected to gauge-in dependent form) and our standard gauge- independent quantities y, yo- The 
equations ( 4.14 - 4.17 ) are generalizations of the formulae ( 2.11 ). 



5 Asymptotics for solutions of Regge- Wheeler equation at 
null infinity 



Let us consider Regge- Wheeler equation (3.40) in null coordinates (u,w 



k 



[du{r^d,y)+d,{r^duy)\ + Ay 



8™ 



-y-0 



(5.1) 



and let us assume (cf. Q) that the solution of (5.1) is in the following asymptotic form (see 
appendix O) : 



E 



(5.2) 



where an are functions of (u, 9, (p) hence they are well defined on J^+. From the assumption that 
an do not depend on v we have for each n: 

d^an = . 

Moreover, denote the u-derivative by dot e.g. 
an — duan • 



Additionally, using (3^) together with the definition of Sa one can express the derivatives of k and 
r 



dyk 



m . 



duk 



dyr = i/c 



dur 



1 



The series (5^) inserted into equation (5T) gives the following formula^ 



E 

n=l 



2(4-n2)a„ 



A flu 



„n+l 



m + 2(n- 1) — ^ H 



o 



= 0. 



(5.3) 



Comparing the coefRcients at the same power of r we obtain recurrence relations for the cocfh- 
cients a„: 



n = 1 ai — free data 

o 

n = 2 2d2 + Aai 



(5.4) 



71 > 2 



2nd„_| 



A + n{n — 1) a„ — 2ma„_i(n — 3)(n + 1) = 



'We assume that the derivatives of the asymptotic terms O ( ^ ) are at least of the same asymptotic order. 
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Let us rewrite equation (5^) using new integer parameter I := {n — 1): 



2(Z + l)d;+2 = - A + ai+i + 2mai{l - 2){l + 2) 



(5.5) 



Remark . The Remark from Subsection 2/7 about NP cons tants remains valid for the Schwarz- 
schild background. More precisely, the right-hand side of (5_^) vanishes on the spherical harmonics 
subspace corresponding to I = 2 and quadrupole part of 04 does not depend on u. 

5.1 "Peeling" for the axial part of Weyl tensor 

We continue our asymptotic considerations based on the assumption ( |5.2| ), in particular the first 
term of the asymptotics gives 



y-0 



The equation (3.38) written in an equivalent form: 

Yd = "^£/(rV)l:a 

may be used to obtain the asymptotic behaviour of ya'. 



(5.6) 



rai + (02 - -ai) + -(03 + mai) + ^(04 + -03) 



1 



(05 — raa-i + 04) + O 



{t y)%v = -2^1 "'"'^i + 2r2 ^ + ^4) + O ( 



or more explicitly, using formula (5_^) we obtain asymptotics of both components of ya'- 



Oi 

r 



From ( 4.17 ) one can show 
r(A+2)Ji?,bcD||A£''V^ 



'2r2 



4m 



(5.7) 
(5.8) 

(5.9) 
(5.10) 



^e,»{r\A"\\Bce^''),b = ^( A + 2)(ry),, + =^(rV),c • (5.11) 



2m , 



From asymp totic s (5.2) for the solutions of the equation ( 3.40 ) we conclude that the -y-component 
of equation (5.11) has the following asymptotic behaviour 



1 / ° s , N 2m , n , / 1 \ 

^(A+2)(ry), + — (rV),.=0(^^j , 

because {ry)^^ ^ (~~) ^ ^ (~) ^^"^ '^^ (5.8). Moreover, the asymptotic condition 



dvX^B = 0(-5-), which is usually fulfilled for the asymptotically flat metric, implies 



4to 
r 



AC\ 



O 
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Hence 



(5.12) 



gives a usual "peeling" because the gauge-dependent term 



-AC 



has the same asymptotics as the full invariant. 
Similarly we may investigate the invariant 



bBd 



TVbdX B 



\\AC 



which equals 
1 



— [ir^yd)tb + [r'^Yb) 



td\ 



(5.13) 



(5.14) 



from ( 4.16 ). Using (|5.2|), (5/7-5^) we obtain asymptotics O(^) for the term ( 5.14 ) if we assume 
that both indices b — d — v. Hence 



2 xu A RC 
r OHv vB\\AC£ 



O 



1 



(5.15) 



Here the result is simpler because SR^AvB is already gauge-independent {-q^y = 0). 

This way we have proved peeling" for axial part of the Weyl tensor corresponding to Newman- 
Penrose scalars and ^'i. In a similar way one can check the asymptotics for the remaining three 
scalars 5*3, ^'4 and ^'2. This is a consequence of the same tensor equations ( 5.11 , ^.13 ) but we 
need to use other components of them. The corresponding invariants have the following asymptotic 
behaviour: 



2)(ry),„ 



2m 



O 



It is relatively easy to verify that the gauge-dependent "correction" terms depending on b have 
the same asymptotic order and finally we obtain 



'^^RabuDWAS 



DA^ab 



XJJ A 

uB\\AC^ 



BC^2 ^ o 



and for the last gauge-independent component 



(5*2 



.ab^ED 



SRa. 



bED 



o 



1 



(5.16) 
(5.17) 

(5.18) 



We would like to stress that asymptotic behaviour of '^2,^3,,^ a given by (5.16-5.18) is not am- 
biguous^. However, in general, and ^'i may have weaker asymptotic behaviour |15[. It is clear 
from the above investigations that axial part obeys "strong peeling" , but this is completely not 
obvious that the same is true for the second degree of freedom described by the polar part. 

^It is shown in Jltj], even for much weaker asymptotic assumptions — so called polyhomogeneous asymptotics 
(i.e. including terms iy(ln r)r~''' where W-polynomial), that the same asymptotics is valid for '^2,'^ 3^'^ i but for 
^0 and "ti, in general, we have terms of order r~^lnr. 
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The results of this section we summarize in the following table: 



Weyl 


C-K-N 


axial invariant 


N-P 


asymptotics 


o ^ 

Wv vB 


2 












o ^ 

W u uB 


9 

' i±CKN 


















+ 2)(ry),„ 










r/3 

— CKN 




+ 2)(ry),„ 




^3 


r-2 




^CKN 




— 2 

r ^y 




*2 


r-3 



Here corresponds to 6R, and in brackets we give the asymptotic results of Christodolou- 
Klainerman-Nicolo (cf. [^). 

6 Conclusions 

In one can find the following statement: "For even waves, it has not yet been possible to 
derive an equation like Regge- Wheeler from the perturbed NP equations" . This question has been 
resolved in [Q but in our opinion not in a satisfactory way (see the discussion at the end of this 
Section). We shall explain in a separate paper why in one can easily formulate Regge- Wheeler 
equation for odd degree of freedom in terms of Newman-Penrose scalars but for polar (even) degree 
of freedom it was difficult to formulate Zerilli equation in terms of Newman-Penrose quantities. 
Although decoupled equations do exist for '^2 and ^'_2 |ltl, even on a Kerr background, they 
are not deformations]^ of the usual wave equations in the asymptotic region. Moreover, the NP 
special null tetrad (chosen in and Q]) is not symmetric with respect to the interchange of null 
coordinates u and v. In other words it has to be chosen in a different way close to future (J^^) 
and past (^~) null infinity. We would like to convince the reader that the Teukolsky equation for 
^pmco _ jg j-^Q^ primary equation describing gravitational waves in asymptotic region (see 
also jl^ for the review of the Teukolsky formalism). The reasons are the following: 

• The Teukolsky equation is not a deformation of a d'Alembert equation in contrast to the 
Regge- Wheeler and Zerilli. 

• The initial data on a slice t — const., instead of position and momenta (x,x), corresponds 
rather to second and third time derivatives of the position (^x, ^x). 

We may think about this equation as an evolution equation for the acceleration x. Usually, when 
we use Fourier transform technique, it is not so important for plane waves which variable (x or 
x) we are using as a canonical position. However, choosing x we exclude from the beginning all 
stationary solutions which are also physically important for the wave operator. 

In our gauge invariant quasilocal formalism we checked the "peeling property" for the various 
components of the linearized Weyl tensor. In this paper we calculated this property only for the 
axial degree of freedom (governed by the Regge- Wheeler equation) and we showed that it is valid 
in its original strong form similar to the case of a flat background (see the table in subsection 
|2.6| ). However, for polar degree of freedom (governed by Zerilli equation |2^, |l^) we may have 
some obstructions. More precisely, because NP scalars are not gauge-invariant one could possibly 
"damage" their asymptotics via gauge transformations. This is related to a more complicated be- 
haviour of the asymptotic solutions of the Zerilli equation and not obvious asymptotics of the gauge 
transformations. Although we showed that for the axial part we could not "damage" asymptotics 

^"By deformation of the wave operator we mean a hyperbolic equation such that when parameter m vanishes it 
becomes the usual wave equation. 
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via gauge transformation, it is not evident that the same is true for the polar degrees of freedom. 
We shall elaborate upon this issue in a separate paper. 

We believe that peeling phenomena for linearized gravitational field is a simpler property for 
the gauge-invariants substituting NP scalars than for NP scalars themselves. 

We also hope that our results can be applied for improving Christodoulou-Klaincrman-Nicolo 
fisll asymptotics on ^+ for nonlinear Cauchy data "sufficiently close" to Schwarzschild. 

It is not easy to relate the results of Q] with our approach. The authors do not give explicit 
formulae for their ^k- However, the deformation of ^'3'° proposed by them differs from ours. In 
our case real and imaginary part of "deformed" do not fulfill the same equation like in 

Moreover, the equations (3.77-3.81) on p. 847 suggest that their invariants ^2, ^2 ^i'*^ quasilocally 
related to 4*4^ and '^q^. We would like to stress that our invariants appear as the "natural" 
positions in the symplectic analysis. The hamiltonian approach enables one to derive the Regge- 
Wheeler equation, resp. the Zerilli equation, as an Euler-Lagrange second order equation for the 
axial, respectively polar, part of the field (see (l^). This suggests that a deformed SR^Pj^ should 
fulfill the Zerilli equation but 3^]/^? fulfills Regge- Wheeler equation. 



A Linearization of Riemann tensor on vacuum background 

To fix the notation and for completeness we present in this appendix the standard linearization 
formulae. Let M be a spacetime with pseudoriemannian metric g^t^. We define a linear perturbation 
of the metric as 

hfiu = Sg^„ gfj,i, - rj^^ . (A.l) 

where by rjf^i, we denote the background metric. The tensor /i^^ is often called perturbation of the 
metric 77^,^. For the inverse metric the perturbation has opposite sign: 

Similarly, we have linearized Christoffel symbols (which are tensors): 

<5r^.A = lrf^iK.;X + Kx;. - /iA.;«) , (A.2) 

where as usual all manipulations are with respect to the background metric, background connection 
etc. We have also linearized Riemann tensor: 

We would like to stress that in general for curved background the linearization of the Riemann 
tensor depends on the position of indices. Let us denote by si^iv '■= /3;ii/ the linearization of 
the Riemann tensor with one upper index, and 

Oi 

Let us compare r^p^^u with the corresponding linearization of the curvature tensor with all indices 
lowered (i.e. SRaup^i = 5 (ganR'^upn))- The relation between 6Raupp and Vaup^ takes the form 

^^avpp ■ — Ravppi^g^ T^av pp{,^i) — ^avpp ^q^^T^ vpp i (-^■^) 

where by Tif upp we have denoted curvature tensor of the background metric. 

We would like to stress that tensor raupp has not all symmetries of the usual curvature tensor. 
This unpleasant property may be verifi ed in the formulae bel ow. 

From definition of r^j/^p and equation ( A.2 ) we can express ( A. 3 ) in terms of h^^, as follows 



^puap c\ ij^pf-pcr ^pv-^ap ^pp-.ua ^up-.pcr ^pa:up ^" ^ua-pp) ■ (-^■^) 
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The basic property of the curvature 



upp 



afip 



impUes 

Let us denote by r^p a hnearization of the Ricci tensor defined as foUows 



' up 



' vap '/ ' pvap 



(A.6) 
(A.7) 
(A.8) 



If we assume that the background metric is Ricci flat (i.e. is a solution of vacuum Einstein 
equations) 

Tlpv = ^, 

then we have the foUowing form for the hnearized Ricci tensor: 
1 



(A.9) 



B Gauge invar iance of y, 



Let us begin with the gauge transformation (B.28) for haB'- 

haB — > haB + ^a,B + ^B.a + 2 EaCs ■ 

r 

and let us denote by Za := /iaS||c£^'^ ^ gauge-dependent axial part of haB- It is easy to check 
that Za transforms as follows: 



Za + {^B\\C£ 



BC\ 



(B.l) 



To compensate the gauge term (C_B||c£^'^),a Ist us analyze the gauge transformation for the axial 
part of XAB = hAB - ^VABhcov'^^ ■ 

XA \\BC^ ^ XA iibc^ +V (^A\\EBC + <iE\\ABc)£ 

Using the commutation relation for the second covariant derivatives 



(B.2) 



^AWEBC = ^AWECB + ^F\\eR^ ABC + £,A\\fR^ EBC ■ 



(B.3) 



where a bc is the curvature of S'^ given by ( 3.18 ), from (B^) and the identity S,^ \\eac^^'" = O7 
we obtain (3.33) 



B ^AC _ B ^AC , f: E AC , ^ <r ^AC 
XA \\BC^ XA \\BC^ + kA\\C E£ +^?A||ce : 

or in equivalent form 

2 B AC 2 B AC , / A I o^^'C ^AC\ 

r XA ||sc£ ^rxA hbc^ + [A +2){^a\\c£ )• 



(B.4) 



Finally the gauge formulae (B.l) and (B.4) imply that ya ~ (A +2)za — {r^A \\bc^ ),a is gauge 
independent. 
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Similarly, we can show gauge invariance of y. We shall derive the relation between y and the 
component rabED^^^ s"'^ of the linearized Riemann tensor which is gauge independent as has been 



shown in Section W. From (A. 7) we obtain 



2rabED — haDtb\\E + hbEta\\D ~ haEtb\\D — hbDtaWE + 

— eb{2haD\\E - '2haE\\D) + — £a(2/l6£;||D - 2/lbD||e) ■ (B.5) 

Moreover from ie^^)ta = ^eaS^^ (cf. (p4|)) we get 

'2rabED£^^ = {haD\\E£^^)tb + {hbE\\D£^^)ta " {haE\\D£^'^)tb " (^b_D| | E£'^'°)ja 

= '2{hbE\\Dl^^^)ta - '2{lT'aE\\D£^^)tb ■ (B-6) 

and finally 

^ED^ab _ Q/i ^ED\ ^ab 

rabED£ £ — ^(llbEWD^ Jta^ 

which implies the demanded relation for y 

y — 2 ^o-bED^ £ • 

This way we get gauge invariance of y from gauge independence of VabED- However, equation 
( 3.36|) also implies this result from the gauge invariance of ya- 



C Compactification of Regge— Wheeler equation near J^^ 

The equation ( |5.l| ) in the coordinates (u,/?, a;'^) can be rewritten^ in the following form: 

2dudp<i> + pdp [kdp{p^)] + A$ + 8mp$ = (C.l) 
where $ := p~^y and k — \ — 2mp. Standard arguments, using domain of dependence considera- 



tions together with conformal covariance of Equation (C.l), show that smooth initial data which 



are compactly supported on some Cauchy hypersurface for the Kruskal-Schwarzschild spacetime 



lead to the solutions of equation (C.l) such that the rescaled $ smoothly extends across J^^ (cf. 



). This means that the assumption (p.2[) is fulfilled for a large class of solutions for the Regge- 



Wheeler equation (5.1). On the other hand, no conditions on initial data which are not compactly 
supported are known, which would guarantee smoothness of solutions across J^+. 

D Axial part of linearized Einstein equations 

The equation ( |3.39| ) contains the component r ^ b\\ac£^^ where by r^^= TS{rAB) we denote 



traceless symmetric part of tensor tab (cf. (3.30)) 

^AB-^ TAB - liVABrcDv'^^ = TS{rAB) , TS{hAB) = XAB ■ 



The general formula ( A. 9 ) takes the following form for tab on Schwarzschild background: 

2rAB = habC"" A^ B + hcoC^ A^ B + h°'B;Aa + B;AF + h"' A;Ba + A-BF (D.l) 
+ hAB-°' a ~ ^AB:^ p — h-AB 

where h := rj^^'^hp^. 

'^^We remind that u := t — r — 2m\n{r — 2m), p := - and are spherical angles. 
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Let us notice that the first two terms in (D^) containing background Riemann C^j^ak are 

o 

proportional to the metric r/AB so they disappear in r^B. To give exphcit expression for tab in 
terms of /i^j^ we need first to examine the terms haB;Ab, hAB-CD and hAB;ab- Using extensively 
the formulae from Section ^ one can show that the following equalities hold: 

hbB-Ad = hhB\\Atd + —^{hABtd^b + "ihfjBWASd — hbftdVAB) + 

^{ihABSdSb + hbfe^EdriAB) + ^{hbdriAB - h-ABSbed) (D.2) 

hAB-CD = hAB\\CD —^"^ {haBWDVAC + hAaWOVBC + haBWCVAD + 

k 

hAaWCVBD + hABtaVCo) + -^{haf£°'£^r]ADriBC + hafs'^e^riBDriAC + 
-hcBVAD - flACVBD - '^hABVCo) (D-3) 

hAB-ab = hABtab + "^—ihABtb^a + hABta£b) + Q^hAB£b£a ^VabflAB (D.4) 

To simplify the analysis we assume that Hab — 0. The final gauge-invariant result is not dependent 
on this assumption but we shall see that formulae are much simpler when we assume this gauge 
condition. Moreover, we may neglect all terms proportional to the metric r]AB because they drop 
out when we pass to the traceless part. From (D.2-D.4) we obtain 

h^B.Aa — h'^BWAfa + '2' — haB\\A^°' 

n- B-AF — -^^—^n-aB\\A£ 

h"' A-Ba — h"" A\\Bta + 2 ^;-/laA| Ise"^ 

n A-BF — —^—^n-aA\\B£ 
hAB;\ ^ 

hAB; p — —£''{'2haB\\A + '2haA\\B) 

h-AB — /i||AS 

where i2 denotes an equality modulo trace: 

tAB—TAB TS{tAB) = TS{tab) ■ 



The above formulae enables one to rewrite (D.l) in a simpler form: 

2rAB - h'^BWAta + h°-A\\BXa - h\\AB ■ (D.5) 

One can easily check the following identity on 5^: 



l-WAB 
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so the last term in ( |D.5| ) drops out when we pass to ^wbc^"^^ ■ Hence 

9° B AC — ( h°- -I- h"- ^ ^DE,a \ \\B ^AC 

AC FB ( la I 7 a DE,a 

= e ri [h B\\A^, + h A\\B^^-VABV h d\\e^. 
Taking into account that 

or simply {r^ri^^)^a — {r^£^^),a = 0, we get 



\\FC 



(D.6) 



2r^rA 



IBC^ 



-AC 



^2^AC^2^FB (/^a^ii^ ^ ^^^^ _ ^^^^^^ h'^ \\pc 



Xa 



(D.7) 



Moreover, the gauge c ondit ion XA \\bc^ = implies that ya = (A +2)ft.aA||c£ and we obtain 



the demanded result (3.39), namely 



The equation (3.38) can be derived in a similar way. Let us start with 



TBd = TadcBrt" + TAdCBri = SRadeBV'"' + SRABEdV T^dS 



(D., 



(D.9) 



where the corresponding components of the linearized Riemann tensor dRf^uXn have the following 
explicit form |]l0| : 

25RadeB — haBXde — haeXd\\B — hdBXae + hcie^a.\\B Q^deWB^a — hae\\B^d) 



H [haBte^d — hdBfe^a + hseta^d ~ h-BetdSa) , 



(D.IO) 



2dR 



ABEd 



•-AdWBE ■ 



hBd\\AE + hBE\\Atd - hAE\\Btd + ^—e'^SdihAaVEB ~ hBaVEA) 



^/fc 

r 



C^hBEWA^d - 2hAE\\B£d + h^dWA^^lBE ~ h ^dWE^" ^1 Ae) 



H — —£°'{hBd%ar]AE - hAdtaVBE + hAatdVEB ~ hBatdVEA) 

Hence from (D.9) we obtain 



(D.ll) 



\/k k 

2rBd — h'^ Btda ~ hdBta"' ^ hdB\\A + 2——{h'^ Bta^d ^ haBtd^"') + ^(hdB — 2£°' SdhaB) 



(hBA^^^-H\\B 



td 



hd'^ta — hd°'xa — -—ha°'ed + h^d\\ 



(D.12) 



Let us notice that the all terms in square brackets (gradient) in the above equality do not contribute 
to axial part rBd\\D£^^ ■ Moreover, if we assume that Hab = (as a gauge condition) we conclude 
that the "second line" in equation (D.12) drops out and the rest gives 



2r 



Bd\\D 



BD 



\fk 1 

= -hdB^^^ADS^^ + Z"|rfa - ZdXa°- + 4 — e°(z<jj„ - Z„|rf) - —Zd 



(D.13) 
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where Za = h^^AWc^^'" (cf- B.l)) and we used the foUowing identity: 



= e 



ED 



\/k k 

h°'BXda — hclBta"" + 2 {h"" BXa^d — haBfdS'^) + 2 — {hdB — ^"'SdhaB) + 2-g/l, 



dB 



impHed by (3.25). Using one more identity 



1 , ° 



which is a straightforward consequence of the commutation relation ( B.3 ) for covariant derivatives 
on 5^, we get 



^BD 



4\/k 1 ° 

[Zatd - Zdja)*" —e^-iZaid " Zd^a) - ^( A + 2)zd 



2rBd\\D£ 
or in equivalent form 

-2r%,|pe^^ = [r\zatb-z,ia)f + r'ya 
The identity 

_ be _ y_ 

^aXd ^dXa ~ ^ad^ ^b^-C — ^ad 2 



similar to (3.36) implies 

-2rVB,||i,e^^-(rV)*''ea6 + rVa 



(D.14) 

(D.15) 
(D.16) 
(D.17) 



which finally gives equation (3.3?) even if we relax the gauge condition hj^g = because both sides 
of ( p.lTl) are gauge independent. 



E Multipoles and traceless tensors 

Let P'' denotes the space of polynomials of degree < fc in M'^. If function / is defined in a 

neighbourhood of a unit sphere S'^, we can denote by Rf its restriction to S'^. 

In Q one can find well-known theorem that R{P'^) is the direct sum X]?=o SH\ where SH^ denotes 

. o 

the space of spherical harmonics of degree /, (g G SH <J=4> Ag = + 1)5). 

Let i be a tensor field in a neighbourhood of in and by TS{Rt) we denote traceless 
symmetric part of Rt on 5^. The following identity holds: 



TS {R{tA^,„A,\B)) =TS {{RtA,...A,)\\B) 



(E.l) 



where "|" denotes covariant derivativ e wit h respect to the flat 3- metric on and "||" stands for 
covariant derivative on S'^. To prove (ET) we can observe the following: 



1 



1 

-VAiBt3A2...Ak + -riA2BtAi3...Ak + 



1 



-VAkBtAtA2...Ak-i3 



and TS ivABtAi...Aj = 0. 

We use spherical coordinates in R'^: 
:= r, (x^), {A = 1, 2), {x^ = 0, x^ 

„ o 

VAB = 7" 1 AB: ms = 1 
^AB 



jAB^ ^33 ^ ^aA ^ q 
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r3 — 1 ^ r A _ 1 



Theorem: 



Proof: 

Let us denote by (x*^) cartesian coordinates in such that 5^ corresponds to the surface: (x^)"^ + 
{x^)^ + (x^)'^ = 1. Transformation rule for transition from spherical to Cartesian coordinates gives 
as following 



P|Ai...A„ =xiAyx'-'^'xPHi. 



ll 12 In 

•'^,Ai'^,A2 ' ' ' ''^ ,AnP^ii---in 1 



and p^i^...in — because degree of the polynomial p is not greater than n — 1. From (E.l) and 
P\Ai...An = we get the result. □ 

In particular for n = 2 we can easily obtain that x^^\\ab is orthogonal to the space SH^ SH^ . 
More precisely, if / G SH° ® SH^ then TS = and 



0= / x^''TS{f^^AB)= I X^^IIAb/ 
Let us consider the following diagram: 



V 



2 '^K yi 11^ V°®V° 

Fl 



where the mappings and the spaces are defined as follows: 

«0l(/,5) = f\\A + £A^9\\B 



ii2{v) = Va\\b + vb\\a - '^ABViic 
i2i(x) = XA^WB 



A AB 

,£ VA\\B 



FKf^ g) = (.9, /) VA = SA^VB XAB = SA^XCB 

- scalars on S*^ 

- covectors on S*^ 

-symmetric traceless tensors on S*^. 

o 

We have denoted by A the laplacian on S . The following equality 

o o 

«io ° «2i o ji2 o «oi = A ( A + 2) 

shows that if we restrict ourselves to the spaces 'Sf := Q [SH° SH^] ( A ( A + 2)V° = V°) 

and V =V^Q [iai{SH^)] (( A + l)V = V ) then the all mappings in the above diagram become 
isomorphisms. 

We define mono- dipole- free scalar as an element of V , mono- dipole- free covector belongs to V 
and any symmetric traceless tensor on 5*^ is mono- dipole- free. 
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